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The  rudiments  of  a  self-consistent  two-phase  thermodynamical  theory  of  intraseasonal  and 
interannual  variability  for  the  tropical  cloud-ocean  mixed  layer  system  are  presented.  In  this  paper  we 
study  some  basic  prop  *ies  of  low-frequency  phenomena  in  the  tropical  coupled  cloud-ocean  mixed 
layer  system  under  Ic  lean  wind  speed  conditions  and  attempt  to  seek  a  unified  thermodynamical 
framework  in  which  tl  ermodynamieal  oscillation  can  be  studied  and  understood.  In  order  to  do  so. 
all  waves  in  the  atmosi  re  and  in  the  ocean  are  initially  filtered  out,  and  the  coupled  system  is  purely 
thermodynamic.  An  ait  ^cean  coupled  model  designed  especially  for  the  low  wind  speed  condition  is 
employed  to  test  the  basic  thermodynamic  feedback  mechanism  between  clouds  and  the  ocean  mixed 
layer.  The  model  has  four  parts:  ilia  shallow-water  system  for  the  large-scale  atmosphere  motion:  (2) 
a  cloud  model.  (?)  a  marine  atmospheric  boundary  layer,  in  which  the  physical  processes  arc 
parameterized  into  the  bulk  formulae  through  a  gcostrophic  drag  coefficient  and  corresponding  heat 
and  moisture  exchange  coefficients;  and  (4)  an  oceanic  mixed  layer  model.  The  coupled  model  is 
solved  analytically  as  an  eigenvalue  problem.  Three  nondimensional  model  parameters  arc  found  to  be 
very  important  in  separating  growing  or  decaying,  oscillatory  or  nonoscillatory  modes:  1 1 1  the  ocean 
surface  stability  index.  < ;  1 2 1  the  surface  wutcr  budget  index,  y,  and  (3)  the  mean  diapvcnal  gradient 
of  the  spiciness  (y)  in  the  entrainment  zone.  5.  The  sign  of  f  divides  ...c  ocean  mixed  layer  into 
shallow  ing  (e  >  <)l  and  entrainment  lr  •  0)  regimes.  The  value  of  y  indicates  the  strength  of  the  w  ater 
budget  versus  the  heat  budget  at  the  ocean  surface.  The  net  freshwater  influx  (e  g.,  precipitation 
exceeding  evaporation)  destabilizes  the  coupled  cloud  and  ocean  mixed  layer.  In  the  entrainment 
regime,  however,  the  parameter  H  becomes  important.  The  unstable  oscillation  only  appears  when  the 
entrainment  zone  is  salinity  dominated.  When  8  exceeds  a  criterion  that  depends  upon  ,  and  v.  the 
oscillatorilv  growing  modes  will  be  generated.  The  model  results  show  that  the  criteria  arc:  8  >  20c  - 
1 1  for  y  —0.1.  8  >  20c  +  1 1.8  for  y  —0.2.  and  >  20c  +  14.6  for  y  -0.5.  The  model  results 
also  demonstrate  that  the  exchanges  of  heat  and  w  ater  across  the  sea  surface  lead  to  both  growing  and 
decaying  modes  of  oscillation  on  tw  o  different  lime  scales,  ow  ing  to  the  stability  of  the  atmosphere. 
For  an  unstable  atmosphere  the  time  scale  is  about  20-30  days.  However,  the  time  scale  is 
approximately  1-3  years  tor  a  stable  atmosphere.  This  work  introduces  the  new  concept  of  two-phase 
thermodynamics  for  tlie  coupled  air-occan  system.  Both  atmosphere  and  ocean  have  two  important 
thermodynamical  variables:  temperature  and  moisture  tor  fractional  cloudiness)  for  the  atmosphere, 
and  temperature  and  salinity  for  the  ocean.  If  salinity  is  neglected  in  the  ocean  model,  no  positive 
feedback  mechanism  will  be  possible  in  the  coupled  air-ocean  system. 


I.  Introduction 

The  ocean  mixed  layer  and  clouds  are  coupled  by  the 
fluxes  of  heal  and  water  mass  at  the  interlace.  The  impor¬ 
tance  of  the  fluxes  of  momentum  and  heat  is  well  recognized 
by  both  meteorologists  and  oceanographers.  However,  the 
water  mass  flux  has  been  given  considerable  attention  only 
in  atmospheric  models,  since  the  latent  heat  release  is  an 
important  source  of  energy  for  the  atmospheric  general 
circulation.  I  he  water  mass  flux  has  been  given  less  atten¬ 
tion  in  ocean  models,  although  it  is  realized  that  evaporation 
and  precipitation  contribute  to  Ihe  surface  buoyancy  flux 
that  influences  the  depth  of  mixing  and  the  thermohaline 
circulation. 

In  a  series  of  papers.  Chit  and  Garwood  [  1989,  1990]  and 
Chu  ci  at.  [19901  presented  positive  and  negative  feedback 
mechanisms  between  clouds  and  the  ocean  mixed  layer. 
First,  clouds  reduce  the  incoming  solar  radiation  at  the 
ocean  surface  by  scattering  and  absorption,  which  cools 
I  relatively)  the  ocean  surface  layer  by  increasing  mixed  layer 
entrainment.  The  cooling  of  the  ocean  mixed  layer  lowers 
Ihe  evaporation  rate,  which  will  diminish  the  clouds.  This  is 
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a  negative  feedback  mechanism.  Second,  precipitation  di¬ 
lutes  the  surface  salinity,  stabilizing  Ihe  upper  ocean  and 
reducing  mixed  layer  deepening.  The  mixed  layer  may  be 
caused  to  shallow  if  ihe  downward  surface  buoyancy  flux  is 
sufficiently  enhanced  by  the  precipitation.  The  reduction  in 
mixed  layer  depth  w  ill  increase  the  sea  surface  temperature 
(SST)  by  concentrating  the  net  radiation  plus  heal  fluxed 
downward  across  the  sea  surface  into  a  thinner  layer.  The 
increase  of  SST  augments  the  surface  evaporation,  which  in 
turn  produces  more  clouds.  This  is  a  positive  feedback 
mechanism.  Figure  I  illustrates  the  feedback  pathways  be¬ 
tween  the  clouds  and  the  ocean  mixed  layer.  1'he  feedback 
mechanism  between  clouds  and  the  ocean  mixed  layer 
depends  largely  on  the  thermodynamics  of  the  ocean  mixed 
layer.  Chu  and  Garwood  1 1990]  and  Chu  el  al.  ]  I990|  predict 
the  excitation  of  oscillatorilv  arid  exponentially  growing 
modes  under  the  strong  surface  wind  ( — 10  m  s  1 )  condition 
w  ith  two  different  time  scales:  (20-30  days.  1-3  years),  that 
largely  depend  on  the  static  stability  of  the  marine  atmo¬ 
spheric  boundary  layer  (MABL).  For  the  unstable  MAUL, 
the  time  scale  is  20-30  days,  a  id  for  the  stable  MABL.  the 
time  scale  is  1-3  years. 

In  the  near-equatorial  region  northeast  of  New  Guinea,  the 
long-term  monthly  mean  wind  speeds  are  of  order  1ms1 
i Sadler  el  aL.  1978).  The  magnitude  of  the  annual  average 
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Fig.  I.  Feedback  paths  between  OPBL  and  clouds. 


vector  wind  stress  for  the  region  4*N-4?S.  1 40- IMF  E  is 
given  by  Wyrtki  and  Meyers  |1976]  as  2.5  x  1 0  •’  N  m;, 
roughly  equivalent  to  a  mean  wind  speed  of  1.2  m  s'1.  It 
remains  to  he  determined  whether  the  feedback  mechanism 
between  clouds  and  the  ocean  mixed  layer  of  Chu  and 
Garwood  [1990]  is  still  important  under  such  weak  wind 
forcing.  In  such  a  case,  the  ocean  surface  buoyancy  flux 
becomes  a  dominating  factor  for  either  generating  or  damp¬ 
ening  of  the  turbulent  kinetic  energy  in  the  upper  ocean, 
depending  upon  the  sign  of  the  buoyancy  flux.  Upward 
buoyancy  flux  causes  the  water  column  to  be  statically 
unstable,  which  generates  turbulence.  On  the  other  hand, 
downward  buoyancy  flux  makes  the  water  column  statically 
stable,  damping  the  turbulence.  In  the  case  of  upward 
buoyancy  flux,  the  mixed  layer  entrains  water  from  below, 
increasing  the  potential  energy.  In  the  case  of  a  strong 
downward  buoyancy  flux  and  weak  winds,  the  turbulence 
kinetic  energy  budget  is  dominated  by  buoyant  damping,  and 
there  is  insufficient  shear  production  to  provide  energy  for 
entrainment.  Hence  the  mixed  layer  is  assumed  to  remain  at 
or  to  be  close  to  a  steady  state,  with  negligible  entrainment 
fluxes  (relative  to  the  surface  fluxes  of  heat  and  water).  It  is 
expected  that  the  effect  of  clouds  on  the  ocean  mixed  layer 
is  different  for  the  stable  and  unstable  regimes. 

Although  our  coupled  model  is  a  one-dimensional  box 
model,  we  are  aware  of  the  importance  of  horizontal  advee- 
tion  and  the  limitations  of  one-dimensional  models.  How¬ 
ever.  the  intent  of  this  work  is  to  develop  a  formalism  to 
examine  thermodynamic  feedback  between  the  two  fluids. 
Because  we  wish  to  concentrate  on  the  thermodynamic 
interaction,  horizontal  advcction  is  either  prescribed  or 
ignored  initially. 


2.  Atmospheric  Model 

The  atmospheric  model  is  a  combination  of  three  parts:  ( 1 1 
a  shallow-water  system  for  the  large-scale  atmospheric  mo¬ 
tion.  (2)  a  cloud  model,  and  (3)  a  marine  atmospheric 
boundary  layer,  in  which  the  physical  processes  arc  param¬ 
eterized  into  the  bulk  formulae  through  a  gcostrophic  drag 
coefficient  and  corresponding  heat  and  moisture  exchange 


coefficients.  A  shallow-water  system  for  the  atmosphere  is 
given  by  [cf.  Gilt.  1980:  Daves.  1985:  I. an  and  Shen.  1988]: 
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where  C\.  and  Va  sire  the  mean  zonal  and  meridional  velocity 
of  the  lower  aimosphcrc.  respectively:  u'a  and  ;y,  are  the 
zonal  and  meridional  perturbation  velocity  of  the  lower 
atmosphere,  respectively:  O',  is  the  perturbation  tropo¬ 
spheric  temperature:  n  is  the  perturbation  fractional  cloud¬ 
iness;  q'  is  the  perturbation  moisture  content,  which  is 
expressed  as  depth  of  liquid  water  per  unit  area  and  related 
to  background  moisture  i/„ :  q,  is  the  saturated  moisture 
content;  I1a  and  represent  the  mean  depth  of  the  lower 
layer  and  the  mean  potential  temperature,  respectively:  S„ 
and  Q'u  represent  the  buoyancy  frequency  and  the  perturba¬ 
tion  heating;  is  the  ocean  surface  evaporation  rale:  P\  is 
the  precipitation  rate:  A,  is  the  total  amount  of  water  vapor 
needed  to  create  a  cloud  over  a  unit  area;  and  JH=  dfidv)  is 
the  latitudinal  gradient  of  the  Coriolis  parameter.  From  the 
mean  distributions  of  temperature  and  moisture  content  in 
the  environmental  air  outside  the  cloud  and  inside  a  deep 
cumulus  cloud  | Kuo.  1965],  we  estimate  that  /t,  —  5  cm. 

The  perturbation  moisture  content  q  and  the  perturbation 
cloud  cover  n  are  computed  for  two  regimes:  convective 
and  nonconvective.  For  the  convective  regime  ( A„  =  I ).  ( \d  I 
becomes  the  prognostic  equation  for  the  fractional  cloudi¬ 
ness: 


w  here  n ’  is  treated  as  the  perturbation  in  fraction  of  clouds 
over  certain  subgrid  areas  in  the  large-scale  atmospheric 
system.  In  order  to  incorporate  the  satellite  data  into  the 
model,  the  subgrid  areas  are  chosen  to  be  1.5°  x  1.5°  (168  km 
x  168  km)  boxes.  The  time  rate  of  change  of  fractional 
cloudiness  is  proportional  to  the  moisture  supply  in  each 
box.  divided  by  the  amount  of  water  vapor  necessary  to 
produce  clouds  in  the  box.  The  main  processes  causing  the 
cloud  dissolution  are  precipitation  and  mixing  with  the 
environmental  air.  The  cloud  evaporation  due  to  mixing  with 
ambient  air  is  a  complicated  problem  and  is  neglected  here 
for  the  sake  of  simplicity.  For  the  nonconvective  regime 
L\a  =  0).  (li/)  becomes  the  moisture  transport  equation: 
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Following  Latt  and  Slwn  [  1988],  the  healing  rate  Q'a  in  the 
convective  regime  is  relaled  to  the  precipitation  rate  by 

Q'o  =  Aut (3) 

where 


v  s  x - ; — 

"aPifipu A  H  a 

which  is  treated  as  a  constant  in  this  research. 

By  linear  regression  of  hourly  rain  amount  and  satellite  IR 
brightness  data  obtained  during  phases  l,  II,  and  III  of  the 
Global  Atlantic  Tropical  Experiment  (GATE).  Albright  ft  til. 
[I985|  suggested  a  linear  relationship  between  precipitation 
rate  averaged  in  the  boxes  and  fractional  cloudiness  /;  of  the 
boxes  by  clouds  with  tops  colder  than  -36°C: 

F,  (ms  ')  =  (0.472  +  8.333/;)  x  |0~7  (4a) 

This  result  verified  Arkin's  |1979|  earlier  analysis  for  the 
GATE  B-scale  array.  From  (4a)  we  obtain  a  linear  relation¬ 
ship  between  perturbation  precipitation  P\  and  perturbation 
fractional  cloudiness 
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dn 

In  order  to  filter  out  all  atmospheric  waves,  we  assume 
b 

—  («„.  0^.  «')  =  6  (5) 

II 

and,  after  using  (3)  and  (4b),  (li/Wh/i  become 
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which  are  basic  equations  for  the  four  unknowns  u'a.  v’a.  0’a. 
and  n‘.  After  solving  (6a)  and  <6 b)  as  algebraic  equations  for 
u'a  and  i>',,  we  compute  the  divergence  of  the  perturbation 
wind  field 


b"L  bti_  Dn 
fi  x  by  D;„  -  /3-y- 


—  - 


ae' 


IK  +  Py  bx 


Wy 


.  *  -T7 


(O',  -  pyK 


ee:, 


be: 


D,n  —  -  py  — 

by  bx 


(7 a) 


If  our  model  is  applied  to  the  equatorial  region,  we  may 
assume  v  —  0  in  (7a),  which  leads  to 


J_ /  ,_  0  3  \  . 
bx  by  ~  Dm  \  l\nTx)*a 


(7b) 


Equations  (6c).  (6</l.  and  (7 h)  are  the  basic  equations  for 
three  variables:  /;'.  and  bu'tl/bx  -  bv'Jby.  Eliminating  the 
last  two  variables  from  (6c).  (6 d).  and  (7b),  we  obtain  a 
diagnostic  equation  for  the  perturbation  fractional  cloudi¬ 
ness  in  the  subgrid  boxes: 


m  \a( 


I 


V  - 


A  a<h  +  U  -  A  „)</,, 


g’  ] 

A«<7<  4"  ( I  —  Aa  )t/(\  J 


A  J>rDmVt 

+ - — 7 —  n'  =  0  (8) 

N-Ul 


where  ih  is  a  second-order  differential  operator  defined  by 
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For  slow  damping  in  the  atmosphere  (i.e..  small  Dr.  Dm).  (he 
fractional  cloudiness  equation  (8)  is  simplified  as 
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In  the  convective  regime  (  =  I),  the  fractional  cloudiness 

equation  becomes 
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where 


(10) 


s  ms(  Tp 

and  if,  is  the  mean  SST.  La  it  and  Slwn  ]  1988]  define  c(  Tr  i  as 
the  SST-dependent  moisture  factor,  w'ith  values  from  0  to  I . 
In  the  convective  regime.  P\  =  in’,  (10)  becomes 

-if <f, in  =  E‘  -  Pi 

With  a  long-time  average,  the  perturbation  model  variables, 
including  fractional  cloudiness  (/;’)  and  ocean  mixed  layer 
temperature  (T,).  should  be  zero,  which  means 


(£*  -  P'r)  dt  =  0 


Xt 


with  a  reasonable  value  of  A/. 


3.  Effects  of  Clolds  on  Ocf:an  Surface  Buoyancy 
Flux 

For  the  case  of  an  ocean  surface  without  ice  (low  and 
middle  latitudes)  the  downward  ocean  surface  buoyancy  flux 
Hi)  has  tw'o  components: 

agF 

B= - 1 - ~p  tj(Pr-E)S  (II) 

P 

Here  a  is  the  seawater  thermal  expansion  coefficient,  p  is  the 
salinity  contraction  coefficient.  p„fl  is  the  characteristic 
seawater  density.  cflH  is  the  seawater  specific  heat  under 
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constant  pressure,  and  F  is  the  upward  heat  flux  plus  net 
radiation  across  the  air-ocean  interface.  The  expressions 
that  follow  are  for  the  net  upward  heat  flux  at  the  ocean 
surface. 


F  -  Rf,-  Rs  +  L,.ph£  +  //,  <  12q> 

for  the  surface  evaporation, 

£  =  PnUuCd.q4Ts)  -  q)/pu  ( 12 6) 

and  for  the  sensible  heat  flux  from  the  ocean  surface. 


Us  *  P'/paU!lCH{Ts  ~  ej  ( 1 2c) 


where  U,,  is  the  gcostrophic  wind.  C„  and  C/,  are  the 
gcostrophic  drag  coefficient  and  heat  exchange  coefficient, 
defined  by  Yamada  1 1976]. 


1 13a) 


(\v» 


Fig.  2.  Dependence  of  Cft.  Ci .  C n  on  the  atmospheric  stability 
parameter  {hJL„). 


only  the  cloud-radiation  and  SST-evaporation  feedback  in 
this  paper.  Therefore 


where  the  von  Karman  constant  is  k  =  0.4.  The  roughness 
parameter  (;„)  is  approximately  1.5  x  1 0 _4  m.  The  ratio 
htJ  Zo  is  0.6  x  to7  for  an  MABL  height.  /»„  =  I  km.  Pr0  is  the 
turbulent  Prandtl  number  lor  neutral  stability,  having  a  value 
of  0.74  according  to  Businperet  at.  [  1971 1.  The  parameters  «. 
6.  and  c  are  experimentally  determined  similarity  functions. 
On  the  basis  of  the  Wangara  data.  Yamada  [  1976]  found  the 
analytical  forms  of  <J.  h.  and  c  as  functions  of  the  instability 
parameter  hJLa,  where  is  the  atmospheric  Obukhov 
length  scale.  It  is  noteworthy  that  the  definition  of  the 
gcostrophic  drag  here  is  different  from  the  quadratic  law 
=  C„U-Y  and  that  the  determination  of  the  analytical 
forms  of  a.  h.  and  c  were  obtained  from  the  Wangara 
atmospheric  boundary  layer  over  land,  rather  than  over 
oceans.  There  have  been  no  comparable  observations  in  the 
MABLas  yet.  Nevertheless,  we  assume  that  these  similarity 
functions  are  valid  for  the  MABL..  Furthermore,  since  the 
computation  method  for  the  moisture  transfer  coefficient  is 
not  well  established,  in  this  study  we  assume  that  CF  =  CH. 
Substitution  of  d.  5.  and  c  into  ( L4o)  and  (146)  leads  to 
apparently  strong  dependence  of  C,;.  C,t.  Q  on  the  atmo¬ 
spheric  stability  parameter  huiL„,  as  shown  in  Figure  2. 
These  parameters  have  much  larger  values  for  the  unstable 
atmosphere  than  for  the  stable  atmosphere,  i.e.. 

C,  -  3. 16x|0  CH,  CE -  2  x  JO hJLa  <  0 

r,-3  x  10  \  C,h  CE-  1.24  x  I0-5  hJLu  >  0 


F' 


bF 

=  —  n' 
bn 


bF 


(14) 


The  effect  of  clouds  on  the  buoyancy  flux  at  the  ocean 
surface  is  twofold:  (I)  decreasing  H  through  the  increase  in 
the  net  heat  loss  at  the  ocean  surface  <F)  by  reducing  the  net 
radiation  at  ocean  surface  (cloud-radiation  feedback),  and  (2) 
increasing  B  due  to  precipitation.  The  variation  of  fractional 
cloudiness  generally  changes  each  term  in  the  right-hand 
side  of  (12a),  i.e..  longwave  back  radiation.  (/?*).  solar 
radiation  (/?,).  evaporation  (£),  and  sensible  heat  flux  <//,). 
Among  them,  only  the  radiation  budget  is  directly  afTected 
by  cloudiness:  the  latent  and  sensible  heat  fluxes  are  indi¬ 
rectly  influenced  by  cloudiness  through  the  change  in  SST. 
Cloud-radiation  feedback  contains  albedo  (reduction  of  solar 
radiation)  and  greenhouse  (emission  of  downward  longwave 
radiation)  effects.  Averaging  the  total  cloud-radiative  forcing 
predicted  by  several  general  circulation  models  | Stephens, 
1989.  Table  2],  we  obtain 


bF  b 

—  =  —  {Rh-R,)~  17  WnT‘  (15) 

bn  bn 


which  means  the  domination  of  albedo  effect  over  green¬ 
house  effect. 

The  variation  in  SST  will  mainly  change  the  longwave 
back  radiation  and  the  sensible  and  latent  heat  fluxes: 
therefore  from  ( 12a).  (126),  and  (12r).  we  can  expect 


The  perturbation  surface  heat  flux  F’  is  computed  by 

bF  bF  bF  bF 

bn  67,  b<i  1  btla 

The  four  terms  in  the  right-hand  side  of  this  equation  indicate 
cloud-radiation  feedback.  SST-evaporation  feedback,  mois¬ 
ture -e\  apteral  ion  feedback,  and  air  temperature-heat  flux 
feedback,  respectively.  Among  them,  the  first  two  feedback 
mechanisms  directly  affect  the  coupled  cloud-ocean  mixed 
layer  system.  In  order  to  simplify  the  system,  wc  consider 


bF 

.  I  I 

67, 


I-I-QS  C  /; 


<  PijR  7,  C  f /  P  I;i  paL,,C  // 


I F  n 


(161 


where  the  Stefan-Boltzmann  constant  is  <rSB  —  5.67  x  10 
W  m~2  K  '4 


-* 


4.  Ocean  Mixed  Laver  Model 

Most  models  that  include  thermodynamic  effects  regard 
the  upper  layer  as  a  well-mixed  turbulent  boundary  layer 
which  exchanges  heal  and  moisture  with  the  atmosphere  and 
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entrains  water  from  below.  The  heat  and  salinity  equations 
take  the  forms 


<‘<7,  F 

hw  —  -  —wr(T}  -  T  h) - +  A  i  ( 17a) 

ht  Pwodpw 

dS 

h*  —  =  -w,(5  -S  *)+(£-  P,)S  +  As  (176) 
dt 

where  /i„,  is  the  mixed  layer  depth;  T and  5_,,  are 
temperature  and  salinity  at  the  base  of  the  mixed  layer, 
respectively:  and  A  ,  and  A  v  are  the  horizontal  advcction  for 
temperature  and  salinity,  respectively.  The  entrainment  ve¬ 
locity  is  w,  and  is  parameterized  by  Chu  and  Garwood  [19X8] 
as 


n>=  A„ 


(C|«i,  -  CzBh) 


fMalT,  -  T-h)  -  01 5  -  5  *)] 


(181 


where  C {  and  C2  are  tuning  coefficients,  and  is  the  water 
surface  friction  velocity,  which  is  computed  by 


^  H'* 


\Pwo/ 


(19) 


Considering  the  low-wind  case  only,  we  set  -  Jems"1. 
uK .  —0.1  cm  s  1 .  The  symbol  A„.  is  a  Heaviside  function  of 
-  CzBlt).  When  (C|«i|.«  —  CzBh)  >  0.  there  is 
sufficient  turbulent  kinetic  energy  to  entrain  and  mix  water 
from  below  and  \„.  I.  which  represents  the  entrainment 
regime.  The  entrainment  rate  is  determined  by  (18)  and  is 
substituted  into  (20)  to  prognosticate  the  mixed  layer  depth 
It- 


(ih„ 

dt 


=  *'<■-  w-h 


(20) 


Here  w  ,,  is  the  mean  vertical  velocity  at  the  mixed  layer 
base. 

When  (Cj«*  —  C2Blt)  <  0.  there  is  insufficient  turbulent 
kinetic  energy  to  entrain  water  from  below,  and  m  ,  is  set  to 
be  zero.  i.e..  A„.  =  0.  This  is  called  the  surface  damping 
regime.  The  mixed  layer  depth  is  calculated  diagnostically 
from  a  balance  of  the  remaining  terms  in  ( 18).  and  it  equals 
the  oceanic  Obukhov  length  scale.  L„ : 


It ..  I- 


C,«J. 

C.B 


(21a) 


The  surface  damping  regime  occurs  when  the  buoyancy  gain 
due  to  excess  precipitation  prevails  over  the  buoyancy  lost 
by  ocean  surface  cooling.  The  mean  mixed  depth  in  the 
surface  damping  regime  is  computed  by 


hw  -  - 


/■'  -3 

C  !«„.» 
C.B 


(216) 


The  temperature  and  salinity  equations  (17 a)  and  (176) 
prescribe  the  balance  of  storage,  entrainment,  and  heat' 
water  mass  flux. 

If  we  use  two  orthogonal  variables  (p„..  ,y)  in  the  TS 
diagram  instead  of  (7,.  S ). 


The  entrainment  velocity  equation  (18)  can  be  simplified  by 

(C,«*„  -  CzBh) 

*»V  =  -A  wPw«  — : -  <18') 

gh(pv  -  p  - 


where  ^  is  “spiciness”  [Munk.  1981].  From  ( 17a).  ( \lb).  and 
(18’),  we  obtain  two  prognostic  equations  for  density  and 
spiciness: 


dp,  ~  CzBhw) 

hw  ~  —  pwn A„ 

dl 


{th» 


P.iyB 

9 


+  p„<((/3j4s  -  aAj)  (23a) 


^X  (C\trw.  -  C.Bhw)  x~x  h  P.  «*> 

hw  ~  —  P«o A„. - - - H - 

ghw  pw-p-k  9 


+  paPA,  +  a  At)  (23  h) 


where 


agF 

D  - - 0</(E  -  P,)S  (24) 

P  i.lTpn 

is  the  spiciness  flux  at  the  ocean  surface,  which  is  orthogonal 
to  the  buoyancy  flux  (11).  The  ratio 


X-X-h  0(5-S_A)  +  a(7,-7-A) 

5  * - =  — - - -  (25) 

Pw  ~  P-h  0(5-5.,,)-  0(7,-  7. „) 

is  the  mean  diapycnal  gradient  of  the  spiciness  [  Veronis. 
1972).  and  is  taken  as  a  model  parameter  in  this  study.  It 
categorizes  the  entrainment  zone:  (I)  6  -  0.  equal  impor¬ 
tance  for  temperature  and  salinity  stratification:  (2)  5  >  0. 
salinity  stratification  dominates:  (3)  S  <  0.  temperature 
stratification  dominates. 


5.  Perturbation  Equations 


The  mean  state  of  the  coupled  cloud-OPBI.  (oceanic- 
planetary  boundary  layer)  system  is  evaluated  from  the 
steady  state  solutions  for  the  equations  (23a).  (23 h).  (18). 
(20).  and  (21  h). 

When  the  coupled  system  is  perturbed  from  its  equilibrium 
state,  the  thermodynamic  feedback  mechanism  between  the 
cumulus  clouds  and  the  oceanic  mixed  layer  either  causes 
the  perturbation  to  grow  (positive  feedback)  or  to  dampen 
(negative  feedback).  The  principal  purpose  here  is  to  study- 
possible  onc-dimensionai  thermodynamic  feedback  mecha¬ 
nisms  between  clouds  and  the  oceanic  mixed  layer,  within 
the  limitations  of  the  simplifying  assumptions.  Hence  the 
energy  exchange  at  the  air-ocean  interface  is  a  primary  focal 
point.  Therefore  we  shall  neglect  initially  the  perturbations 
of  those  variables  not  directly  related  to  the  exchange  at  the 
air-ocean  interface,  such  as  the  perturbations  of  the  horizon¬ 
tal  advection  A'T  and  A's.  the  perturbations  of  the  atmo¬ 
spheric  variables  at  the  top  of  the  MABL,  q'  and  O',,,  and  the 
perturbations  of  the  oceanic  variables  at  the  base  of  the 
OPBL,  hT,,.  71, ,.  and  S'.h.  From  the  basic  equations  of  the 
coupled  system  (23a).  (23 b).  (18').  and  (20)  or  (21a).  the 
equations  for  p'v.  x  arc 


hpi  _  Pud  1  +  C: Aw) 
dt  ()RW 


Pw  =  PhO<05  -  «7)  *  =  p„<t(05  +  a 7)  (22) 


(26a) 
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i>K  Cip„<t/?  hi  CipKil 
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^  </Ap„  /i„  Ap„ 

where 


Ap„.  =  p„.  -  P-A 


TABLE  I-  Values  of  The  Model  Parameters 


Parameter 

Value 

IK, 

1 .29  kg  m  '3 

h< 

5  cm 

u„ 

3  ms"1 

Pu 

1035  kg  m  ~3 

T, 

28  C 

.V 

35  g  kg'1 

/T„ 

50  m 

a 

0.2  x  10  ‘  K  1 

p 

0.8  x  I0'} 

C, 

1.0 

c. 

0.2 

is  the  mean  density  jump  at  the  base  of  the  mixed  layer,  and 
B’  and  /)’  are  perturbations  of  surface  buoyancy  and  spici¬ 
ness  fluxes,  which  are  computed  by 

T ;  -  P&P,  -  E)S •  (276) 


If  the  low  surface  wind  condition  is  considered  (especially 
for  the  western  Pacific  warm  pool  region),  it  is  reasonable  to 
assume  that  the  upwelling  velocity  (w_,,)  is  also  a  negligible 
quantity.  We  may  drop  the  last  term  on  the  right-hand  side  of 
(26c)  for  the  first  approximation. 

The  definitions  of  surface  buoyancy  and  spiciness  fluxes 
(ID  and  (24)  lead  to 


ag  <)!■  _  bP, 

- —  4 -figS  — 

PwOFp* 


an  bF 

Pn&pv  bn 


(I  -  A)  (28 a) 
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ag  bF  _  bP,  ag  bF 

-  p<jS  —  - - - —  —  f  I  +  A) 


P.flfp,.  bn  bn  p,r(K>„.  bn 

bB  ag  bF  bE 

—  = - - - pgS  — 

bT%  Pn<)C/>w  bT %  b /  , 

bD  ag  bF  _  bE 

—  = - +  pgS  - — 

bl  -x  PuflCpu'  b  T$  b  7\ 


(28  h) 

(28c) 

(28J) 


where 


A  = - = - 

abF/bn 

From  (41.  (15),  and  Table  I.  the  parameter  A  is  estimated  as 
A  -  0.034. 


6.  Time  Scales 

Chu  el  al.  [1990]  found  three  lime  scales  in  the  thermody¬ 
namically  coupled  oceanic  and  atmospheric  boundary  lay¬ 
ers. 


6.1.  Fractional  Cloudiness-Damping  Time  Scale 


1  bPr 
h,  bn 


(29) 


Using  (4)  and  taking  lic  -  5  cm.  we  have  r„  -  0.6  day.  Since 
n  is  interpreted  as  fractional  cloudiness  averaged  over  sub¬ 
grid  areas,  the  spatial  scale  for  the  fractional  cloudiness- 
damping  time  scale  is  around  168  km. 


6.2.  Ocean  Mixed  Layer  Cooling  Time  Scale 


I  bF 
P  *,C  „  bT\ 


(30a) 


On  the  basis  of  the  assumption  that  ba  and  q  in  (126)  and 
(12c)  are  determined  by  the  large-scale  atmospheric  motion 
only,  from  the  surface  heat  balance  equation  (I7«)  we  have 


Tr* 


L;xp  cE 

CpvR  T'  C/i 


P<A  ;><:  I-  H 


C„ 


(306) 


6.3.  Cloud-SST  Coupling  Time  Seale 


,-t  _ 

n.T 


P»-cPjTu.  bn  bT ; 


bF  bE 


This  lime  scale  r„  T.  is  the  time  needed  for  the  process  that 
the  deviation  of  fractional  cloudiness.  n‘ .  due  to  the  change 
of  ocean  surface  evaporation  caused  by  SST  variation,  feeds 
back  to  SST  through  its  effects  on  the  surface  heat  flux. 
Utilization  of  (15).  (16).  and  (27)  leads  to 


f  T„U,jLq,[biRh  -  /?,)/an]l 
[  P  pi-t  ffvlt  i*>6  rR  i>T  i 


(31) 


Among  these  three  lime  scales.  t„  is  the  shortest.  The 
other  two.  r,  and  r„  , .  depend  on  the  parameters  C..  Cq. 
and  Ct.  which  arc  functions  of  the  atmospheric  stability. 
Figure  3  shows  the  dependence  of  rT  and  r„T  on  the 
atmospheric  stability.  For  the  unstable  atmosphere  ( hJLa  < 
0).  7„r  -  3-6  days,  and  rr  —  20-30  days.  For  the  stable 
atmosphere  OiJL„  >  0),  r„  T  ~  0.3-1  years,  and  t-,  -  1-3 


years. 

The  mean  seasonal  outgoing  longwave  radiation  (OLRi 
consists  of  contributions  from  higher-frequency  components 
of  convection.  Lau  and  Chun  [1988]  found  two  bands  (1—5. 
40-50  day )  of  oscillations  in  OLR.  Large  variability  in  the  I- 
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Fig.  3.  Dependence  of  -/  and  r„  ,  on  the  atmospheric  stability 
parameter  </i(, //.„). 


7.2.  Surface  Water  Budget  Index 

The  relative  importance  of  surface  freshwater  influx  in  the 
surface  buoyancy  flux  can  be  measured  by  the  following 
parameter: 

Pf)S(Pr  -  E) 

g—  <33) 

7.3.  Mean  Diapycnal  Gradient  of  Spiciness  at  Mixed 
Layer  Base 

This  parameter  appears  only  in  the  entrainment  regime. 
The  definition  of  this  parameter  (5)  is  given  by  (25).  It  is  a 
tracer  conserved  by  isentropic  motions  and  an  indicator  of 
diffusive  stability.  It  should  be  useful  for  the  combined 
description  of  heat  and  salinity  transfer  in  the  entrainment 
zone. 

8.  Non  dimensional  Perturbation  Equations 


to  5-day  band  is  confined  with  the  cloudy  region  distribution 
(e.g..  equatorial  western  Pacific).  The  coincidence  of  these 
two  bands  with  - , ,  and  rr  for  the  unstable  atmosphere 
tintense  convection)  indicates  that  rT  (20-  to  30-day  time 
scale)  is  one  of  the  intraseasonal  variations. 

The  spectrum  of  average  SST  anomalies  near  the  upper 
Peru  Coast  (4'-I26Sj  exhibits  a  major  peak  at  36.6-42.7 
months,  and  a  second  peak  at  10  months  [Rasntusson  and 
Carpenter ,  1982].  The  MABL  near  the  Peru  Coast  is  usually 
stable.  The  concurrence  of  these  two  peeks  with  rn  T  <  100— 
300  days)  and  rT  (1-3  years)  for  the  stable  atmosphere  may 
imply  that  our  theoretical  model  is  applicable  to  the  low- 
frcquency  modes. 


By  using  (28)  and  the  definitions  for  the  time  scales  (rn  /. 
tj).  t30a).  (31).  and  for  the  model  parameters  y.  the  pertur¬ 
bations  of  surface  buoyancy  and  spiciness  fluxes  are 

B'  ~  ~ —  17t  '  +  V/0  “  A)](*'  “  Pie) 


By 


2p»oPS 


p(*'  +  pi)  (34  a) 


9>i  H 
2p»« 


+  Vt<1  -  A)]0r'  “  Pi) 


By 

2pwap$(* 
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(34  b) 


7.  Model  Parameters 

Three  nondimensional  parameters  are  found  important  in 
the  coupled  system. 


7. 1 .  Surface  Stability  Index 

The  ocean  mixed  layer  exhibits  different  behavior  ten- 
training  or  shallowing),  depending  on  the  sign  of  the  param¬ 
eter: 


7n.tiSfiH  —  C  |//;|  JC2) 

Prfis 


(32«) 


When  e  >  0.  the  ocean  mixed  layer  is  in  the  surface  damping 
regime.  When  e  <  0.  the  ocean  mixed  layer  is  in  the 
entrainment  regime.  In  this  study,  we  deal  with  the  very'  low 
surface  wind  condition,  which  leads  to 


C;/)» 


From  (33)  the  ratio  of  the  mean  surface  spiciness  flux  to 
mean  surface  buoyancy  flux  is 

D 

—  =  1  -  2y  (35) 

B 

Time  t  and  the  perturbations  are  nondimensionalized  by 

/  _  _  yhl.  _\ 

{1.  p„..  /1*,)  —  |  ” 7 1 >  P*oP**,  P hH.J 

(36) 

Assuming  the  upwelling  velocity  to  be  negligible,  the  nondi¬ 
mensional  form  of  the  perturbation  equations  (26 a).  (26 b). 
and  (26c)  are 

3p*.  I  +  C2A*.  1  +  C2A* 

—x  - - - - (I  +  ma  +  )p».  + - - - (I  +  ma  )x 

+  (I  +  C2A*.)m/i*  (37a) 


Therefore  the  surface  damping  factor  is  approximated  by 

T/i.r# 


P<jftnS 


(32  b) 


9X  1  , 

—  =  -{ 1  +  ma  -  C2A„.S(  1  +  ma  .  )]p„. 
dt  2 

— 5  f  1  +•  ma  +  -  C2A*.S(  1  +  mu  _  )]* 


which  means  that  the  surface  buoyancy  flux  is  dominant 
over  the  surface  wind  forcing. 


-  (I  -  2y  -  C2Ah.S  1/71/1*  (37(7) 
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Fig.  4.  Separation  of  modes  in  the  Ic,  y)  plane  for  the  surface 
damping  regime  (r  -•  0). 


M  dfilk  C2 
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where 


a,  =  I  —  A  +  ey  «  .  =  1  —  A  —  ey 


~t  Ap„, 

m= -  ix  = - -  (38) 

TnJ  Ph«PS 

Only  three  parameters  in  the  basic  system  (37a).  <37/>).  and 
(37c)  are  allowed  to  vary:  e.  y.  and  (5.  The  others  are  fixed: 
A  -  0.034.  m  ~  3-5,  and  p. - 0.01. 

9.  Solutions 

The  system  (37a).  (37/>).  and  (37r)  contains  both  surface 
damping  (A„  =  0)  and  entrainment  (A„,  =  I)  regimes.  The 
general  solutions  of  the  basic  system  (37a).  < 37/>».  and  (37c) 
have  the  form 

<Pw.  X.  Kl  ~  /V  (39) 

where  £>v.  !),„{  j—  1. 2.  *  •  •  )  are  the  integral  constants, 
and  <r,  are  eigenvalues,  which  are  roots  of  an  algebraic 
equation.  It  is  noteworthy  that  for  the  entrainment  regime 
(A„  =  I).  this  algebraic  equation  is  of  the  third  order, 
whereas  for  the  damping  regime  (  \u.  =  0).  it  is  of  the  second 
order.  This  order  change  is  attributable  to  the  change  from 
prognostic  to  diagnostic  form  for  the  respective  mixed  layer 
depth  equations  (37 c).  Therefore  in  the  solution  (39).  there 
arc  two  terms  for  the  surface  damping  regime  and  three 
terms  for  the  entrainment  regime. 

10.  Instability  and  Oscillation  Criteria 

The  instability  criterion  for  the  thermodynamically  cou¬ 
pled  air-ocean  system  is 

<rr  =  Re  (<r)  <  0  decaying 

it,  =  Re  (o-)  =  0  neutral  (40) 

a,  =  Re  (it)  >  0  growing 

The  oscillation  criterion  for  the  coupled  system  is 


Fig.  5.  Eigenvalue  <r  in  the  !<-.  y)  plane  for  the  surface  damping  regime  ( e  >  0):  (u)  trr;  (b)  | at . 
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<r,-s  Im  (<r)  =  0  nonoscillatory 

(41 ) 

<r:  =  1m  (<r) 5*0  oscillatory 

We  should  keep  in  mind  that  the  linear  theory  is  valid  only 
for  infinitesimally  small  perturbations,  i.e..  j(  )'|  «:  |(~) .  As 
the  perturbations  become  sufficiently  large,  this  requirement 
is  no  longer  satisfied:  therefore  the  solutions  are  no  longer 
valid.  We  must  then  consider  the  nonlinear  system. 

10. 1 .  Shallowing  Regime  (e  >  0) 

In  this  regime  (A...  =  0).  substitution  of  (39)  into  (37 a), 
Gib),  and  (37c)  leads  to  a  quadratic  equation  of  <r: 

[cr  (1  +  +  —  [I  +  ma+  -  ( I  t  ma  .  )(1  -  2y)]  j 

I 

+  ~  [1  +  ma~  +-  (I  +  mc+)(l  -  2y)]ll  ma  ( =  0  (42) 


When /((c)  >  y  >  /,(e).  the  positive  and  negative  feedback 
mechanisms  arc  comparable.  Because  of  phase  differences, 
the  corresponding  modes  are  oscillatory.  If  y  is  sufficiently 
large,  i.e..  y  > />(c).  the  positive  feedback  mechanism  caused 
by  the  net  surface  freshwater  influx  prevails,  and  the  system  is 
oscillatory  and  unstable.  The  growth  ranges  from 
0  —  I  r  1  (Figure  5 a),  and  the  period  of  the  oscillation  is  the 
order  of  27tt7  (Figure  5b).  However,  for  y  < /.< c),  the  negative 
feedback  mechanism  dominates,  and  the  initial  perturbations 
will  be  damped.  Furthermore.  (37 c)  shows  that  in  the  entrain¬ 
ment  regime  (A„  =  1).  if  the  density  jump  at  the  mixed  layer 
base.  Ap„  ,  is  infinitesimally  small,  the  local  time  rate  of  change 
of  the  mixed  layer  depth  will  be  negligible.  Equation  (37c)  is  the 
same  as  for  the  surface  damping  regime.  Therefore  when  the 
mixed  layer  depth  is  predicted  to  increase  by  (37<  )  for  the 
surface  damping  regime  (  \H  0).  the  regime  switches  to  the 
entrainment  regime,  with  negligible  density  and  spiciness 
jumps  assumed  at  the  base  of  the  mixed  layer. 


It  is  notable  in  (42)  that  the  parameter  S  does  not  appear  in 
the  surface  damping  regime.  There  are  two  roots  of  cat  each 
pair  of  given  values  of  model  parameters  (e,  y).  Between  the 
two  roots,  the  first  one  always  has  a  negative  real  part 
through  the  whole  r-y  plane,  representing  damping  modes, 
w  hich  arc  not  of  interest  and  are  neglected.  The  second  root 
has  a  positive  real  pan  in  some  domains  of  the  t-y  plane.  The 
properties  of  this  root  indicate  different  modes  according  to 
its  real  (a>)  and  imaginary  (tr;)  pans 


For  the  low  wind  speed  ease,  the  inequality,  e  <  0.  is 
approximately  equivalent  to  8  <  0.  The  parameter  y  is 
negative  (positive)  when  the  mean  precipitation  rate  (Pr) 
exceeds  (lags)  the  mean  surface  evaporation  rate  (E).  The 
basic  system  l37u).  Gib),  and  (37c)  of  this  regime  is  a  set  of 
three  first-order  ordinary  equations.  Substitution  of  (39)  into 
(37«).  (37/j).  and  (37c)  leads  to  (A.  =  l)a  third-order 


IXr  >  0 

07  =  0  =>  nonoscillatorily  growing 

algebraic  equation  for  the  eigenvalue  <r: 

a>  >  0 

<r,^0  =>  oscillatorily  growing 

(43) 

An  —  cr  A 12  A  u 

<7,  <  0 

rr,sQ  =>  oscillatorily  decaying 

A 12  A22  “  <7  A  23. 

A 

0 

<r,  =  0  ^  nonoscillatorily  decaying 

A31  A32  A  33  —  a 

-  0 


(46) 


Figure  4  shows  the  separation  of  different  modes  in  the  (c-y) 
plane.  The  three  lines  (if  small  fluctuations  are  ignored) 
divide  the  U.  y)  plane  into  four  modes: 


y  >/|(e)  =>  nonoscillatorily  grow  ing 
/i(e)  >  y  >fyr  )  oscillatorily  growing 
/tU)  >  y  >/df  I  oscillatorily  decaying 


(44) 


where 


An  -  — - 


A 1 2  — 


(1  —  ma.) 


(l-i-  mu  .) 


y  </)(«)  nonoscillatorily  decaying 

where 

<*•  2  +  I5e  file)  =  2  +  6e 

(45) 

fye)  =  2  -  20r(e  <  0.05)  /,(*-)  =  I(e  >  0.05) 

The  model  parameter,  y.  represents  the  relative  importance 
of  the  net  surface  freshwater  influx  in  the  mean  buoyancy 
flux  B  and  indicates  the  strength  of  positive  feedback  mech¬ 
anism  between  the  coupled  system.  On  the  other  hand,  the 
surface  damping  factor  (e)  denotes  the  negative  feedback 
mechanism.  When  y  >  the  positive  feedback  mecha¬ 
nism  dominates  the  negative  feedback  mechanism,  and  the 
prognostic  variables  (p„ .  g)  increase  exponentially  with  the 
growth  rate  ar  --  3-5  (in  units  of  rf ').  The  time  during  which 
the  perturbations  of  density  and  of  the  spiciness  <>) 
double  their  strengths  is 


A  id  —  ( I  +  ( 2.)w 

1 

Aji  =—  [I  +-  ma  —  CjSd  +  ma-)] 

(47) 

1 

A 22  =  — —  [  I  +  ma  -  —  CiSIl  +  ma.)] 


Ajs  =  —  (1  —  2y  —  C2S)m 
C2e 

A31  =  — —  (I  +  ma-) 
C2e 

A  i2  - - ( I  +  ma  ) 

2m 


!n  2 

7  double  —  0.2t  f 

<*r 


C-yum 
Ajj  =  — - 

M 
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Fig.  6.  Separation  of  modes  in  the  (r,  R)  plane  for  the  entrainment 
regime  (c  <  0)  for  y  =  -0.1. 


There  are  three  roots  of  <r  for  each  set  of  model  parameters 
(c.  y.  8).  Among  the  three  roots,  the  first  always  has  a 
negative  real  part  and  represents  damping  mode  that  is  not  of 
interest  here.  The  other  two  roots,  in  certain  domains,  are 
complex  conjugates,  representing  the  same  oscillatory 
modes.  In  other  domains,  they  arc  both  real,  representing 
the  nonoscillatory  modes.  The  eigenvalue  with  a  larger  real 
pan  will  be  discussed  in  this  section.  The  distributions  of  tr, 
and  <r,  control  the  properties  of  the  modes,  as  shown  in  (43). 
The  solutions  show-  that  only  damping  modes  are  available 


O 


Fig.  7.  Separation  of  modes  in  the  I  r.  a)  plane  for  the  entrainment 
regime  tr  <  0)  for  y  =  —0.2. 


when  the  surface  evaporation  rate  E  exceeds  the  precipita¬ 
tion  rate  Pr.  In  the  low  mean  wind  speed  region,  the 
entrainment  process  is  mainly  driven  by  upward  buoyancy 
flux  (i.c..  B  <  0).  From  the  definition  of  y  (33),  the  bigger 
absolute  value  of  y  implies  the  weaker  entrainment  process. 
For  a  not  very  weak  entrainment  process  (i.e..  |yj  is  not  too 
large),  the  solution  greatly  depends  on  the  mean  diapycnal 
gradient  of  spiciness  in  the  entrainment  zone.  6.  We  compute 
the  eigenvalues  in  the  (e,  S)  plane  for  several  negative  values 
of  y  -0.1.  -0.2.  and  —0.5.  as  shown  in  Figures  6-8.  In  the 
entrainment  regime  (c  <  0).  the  net  freshwater  influx  makes 
y  <  0.  In  such  a  case  (f<0.  y<0).  the  parameter  S  will  play 
a  key  role  in  determining  the  stability  of  the  coupled  system. 
When  8  exceeds  a  criterion  that  varies  with  f  and  y.  the 
oscillatory  growing  modes  will  be  generated.  The  model 
results  show'  the  criteria  to  be  8  >  20s  +  1 1  for  y  =  —0. 1 
(Figure  6),  5  >  20c  f  1 1.8  for  y  =  —0.2  (Figure  7).  and  S  > 
20c  +  14.6  for  y  =  -0.5  (Figure  S).  The  shrinkage  of  the 
unstable  oscillation  area  with  the  decrease  in  ytfrom  -0. 1  to 
-0.5)  implies  that  (I)  the  more  the  net  freshwater  influx,  or 
(2 1  the  weaker  the  entrainment  process,  the  stronger  the 
salinity  stratification  domination  should  be  in  the  entrain¬ 
ment  zone. 

The  distributions  of  the  real  and  imaginary  parts  of  the 
eigenvalue  <rr.  ioy|  for  different  values  of  y  are  shown  in 
Figure  9  (y  =  —0.1),  Figure  10  (y  =  —0.2),  and  Figure  II 
(y  -  -0.5).  If  our  interests  are  only  in  the  unstable  modes, 
representing  the  enhancement  of  initial  perturbations,  these 
figures  indicate  the  following: 

I.  In  the  unstable  region,  the  decrease  of  parameter  y 
(from  —0.1  to  —0.5)  decreases  the  growth  rate  tr,  and 
increases  the  frequency  |trj  (1.5  to  3.5  7-j  ').  This  indicates 
that  the  ability  of  the  net  freshwater  influx  (represented  by 
the  negative  value  of  y)  to  amplify  the  initial  perturbation  is 
reduced  when  the  entrainment  process  weakens. 


(Jjci/liilnrily  Gromlnx 


Fig  8.  Separation  of  modes  in  the  (c.  R)  plane  for  the  entrainment 
regime  (c  <  0)  for  y  =  —0.5. 
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Fig.  9.  Eigenvalue  a  in  the  (t.  5)  plane  for  the  entrainment  regime  If  <  0)  and  lor  y  =  -0.1:  (<i)  ar.  and  (b)  |oj[. 


2.  The  unstable  oscillation  only  appears  as  5  >  0,  which 
means  that  the  net  surface  freshwater  influx  generates  the 
unstable  modes  only  in  the  salinity  stratification-dominated 
entrainment  zone. 

3.  When  the  entrainment  zone  is  dominated  by  temper¬ 
ature  stratification,  the  net  surface  freshwater  influx  can  only 
dampen  the  initial  perturbation. 

1 1 .  Conclusions 

The  feedback  between  the  cloud  and  OPBL  in  the  coupled 
MABL  and  hydrodynamicallv  stable  OPBL  system  is  inves¬ 
tigated  by  a  simple  one-dimensional  coupled  model.  The 


time  scales  largely  depend  on  the  stability  of  the  MABL.  For 
the  stable  atmosphere,  the  two  time  scales  are  quite  long: 
rn.r  ~  100-300  days,  and  tt  —  1-3  yr.  For  the  unstable 
atmosphere,  however,  the  two  lime  scales  are  much  shorter: 
~n.r  ~  3-6  days,  and  tt  ~  20-30  days.  In  the  western  Pacific 
warm  pool  regions,  the  MABL  is  usually  unstable.  Therefore 
this  theory  may  provide  some  explanation  for  the  two  lime 
scales  (3-6  and  20-30  days)  of  intense  atmospheric  convec¬ 
tion  above  the  western  Pacific  warm  pool. 

The  mean  freshwater  influx  at  the  ocean  surface  due  to  the 
excess  precipitation  over  evaporation  amplifies  the  initial 
perturbations. 


Fig.  10.  Eigenvalue  cr  in  the  (r,  S)  plane  for  the  entrainment  regime  U  <  0)  and  for  y  -  -0.2:  («)  <r,.  and  (h)  |er;|. 
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Fig.  II.  Eigenvalue  >r  in  the  \f.  tit  plane  for  the  entrainment  regime  tr  <  0)  anti  for  y  -0.5:  <«>  <r,.  and  <J»  |<r/. 


In  the  enirainment  regime,  the  mean  diapycnal  gradient  in 
the  enirainment  zone  <5  plays  an  important  role  in  determin¬ 
ing  the  instability  of  the  coupled  system.  When  8  exceeds  a 
criterion  that  varies  with  f  and  y.  the  oscillatorily  growing 
modes  will  be  generated.  The  model  results  show  that  the 
criteria  arc  §  >  20f  +  1 1  for  y  —  -0.1.  6  >  2l>r  -  1 1.8  for 
y  =  -0.2.  and  8  >  20f  -  14.<S  for  y  -  -0.5. 

Furthermore,  the  dynamics  of  the  ocean  mixed  layer  is 
highly  nonlinear.  There  is  no  smooth  transition  between  the 
entrainment  and  the  surface  damping  regimes.  The  use  of  the 
linear  theory  enables  our  initial  search  for  feedback  mecha¬ 
nisms.  The  nonlinear  effects  should  be  take  into  account  in 
realistic  applications. 
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